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Abstract 

Climate change is a reality of today. Paleoclimatic proxies and climate predictions based on 
coupled atmosphere-ocean general circulation models provide us with temperature data. Using 
Detrended Fluctuation Analysis, we are investigating the statistical connection between the climate 
types of the present and these local temperatures. We are relating this issue to some well-known 
historic climate shifts. Our main result is that the temperature fluctuations with or without 
a temperature scale attached to them, can be used to classify climates in the absence of other 
indicators such as pan evaporation and precipitation. 
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I. INTRODUCTION 



Our climate is changing rapidly. We see the effects of this climate shift in every aspect of 



to changing migration 
8fl. While the scientific 



nature from retreating glaciers [l, and breaking ice shelves {3], 0, 
patterns of animals g, |7| and even the extinction of certain species 
evidence for a serious climate shift is growing, there are still some groups trying to cloud the 
public opinion with arguments about past climate changes having little or no effect on our 
society. The first rudimentary thermometer was invented by Galileo Galilei in 1593. Daniel 
Gabriel Fahrenheit invented the first mercury thermometer and made it possible for scientists 
to measure temperatures consistently. Our knowledge about the climates before the time of 
Galileo and Fahrenheit usually comes from written (or sometimes painted) personal accounts 
or what kind of wine was produced in that region of the Earth, rather than from irrefutable 
scientific evidence. Although many such accounts are very useful in determining the climate 
of the "Old World", we have settled many different parts of the Earth, and hence more 
information is needed about paleoclimates to make more sense of the changes happening 
today. 

Our scientific efforts are normally based on the analysis of temperature records obtained 



from various sources like tree rings [9|, ice core data [10|], ocean or lake sediments [ll|, Il2j . 
or coral reefs (l^ . Since these proxies provide us mostly with a temperature record, this 
only gives us limited knowledge of the climate. For example, if we compare Ospina Perez, 
Colombia to Las Vegas, USA, we can immediately say that the climate of these two cities 
should be vastly different. Ospina Perez, Colombia has a tropical monsoon climate receiving 
about 15 times more rain than Las Vegas, USA which is located in a dry/hot low latitude 
desert. However, if we look at the yearly average temperatures of these two cities, we find 
the averages to be around 19.1 o C. Therefore, in addition to temperature records, the 
climatologists are also using average precipitation and evaporation rates to determine the 
climate of a certain region. 

Determining the weather is actually a rather simple issue. A cold day is usually followed 
by a cold day, and a warm day is usually followed by a warm day. On a larger scale, a colder 
week is usually followed by a warmer week which corresponds to the average duration of the 
general weather regimes. But as the longer timescales are governed by different processes 



like circulation patterns and global warming, defining long-term correlations becomes more 
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difficult. Recent 
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V, scaling arguments have been used to ana 



In order to separate the trends and the correlations we need to eliminate the trends in 



our temperature data. Several methods are used effective 
3 3 44 



analysis (R/S) 
analysis (DFA) 
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451 ] . wavelet techniques (WT) 
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or this purpose: rescaled range 



461 ] and detrended fluctuation 



47|. 



Climate change is a major concern for everyone, and it is a focus of atmospheric research. 



The changes in properties of the cloud cover influence the radiative feedback 



and become one of the important players that effect the climate in general 4l|. Considerable 



rom the clouds 



effort has 
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3een made to study the time scaling of fluctuations of the cloud cover using DFA 
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Analysis of the temperature fluctuations over a period of decades in different places of 
the globe has already shown the effectiveness of the application of DFA to characterize the 
persistence of weather and climate regimes. DFA and WT have been applied to study tem- 
perature and precipitation correlations in different climatic zones on the continents and also 
in the sea-surface temperature of the oceans. The recent results show that the temperatures 
are long range power law correlated. The long-term persistence of the temperatures can be 
characterized by an auto-correlation function, C(n), of temperature variations where n is 
the time between the observations. The auto-correlation function decays as C(n) ~ n -7 . 
Though there is some disagreement on the value of 7, the fact that the persistence of the 
temperatures can be characterized by this auto-correlation function is firmly established. 



Different groups have used R/S, DFA and WT ana 
same value, 7 ~ 0.7 for the continental stations 15 . 



ysis to s 
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0.4 for island stations 
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18| and on the oceans 
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low that_7 has roughly the 



M lM [23, and roughly 



181 ] . These methods have also 



been applied to the temperature predictions of coupled atmosphere-ocean general circulation 
models [19, 30, [4^, [4^] but there is disagreement on the actual value of the exponent 7. On 
one side it is argued that the exponent does not change with the distance from the oceans 

n 

|15l ] and is roughly 7 ~ 0.7. On the other side it is said that the scaling exponent is roughly 
'.. over the oceans, roughly 0.5 over the inner continents and about 0.65 in transition regions 

Previous work in this area has also shown that there is a slight variation in the sea. 



exponent between the low-elevation, mountain, continental and maritime stations 27 
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43|. Most of these variations are within the standard deviation of 7 and even though some 



of the previous work shows a correlation with location and elevation |27l . I28j | it is difficult to 
establish a general relationship between the statistical nature of the temperature fluctuations 
and elevation. 

In our previous work we have learned that Detrended Fluctuation Analysis of temperature 
data can be used to suggest a relationship between the temperature fluctuations and climate 



without any help from precipitation and evaporation data [42], |43( . The present work tries to 
extend the same method to show that DFA can be used to discover historic climate changes 
for which no scientific climate records exit. The natural extension of this work will be to 
quantify present climate changes. 

II. METHOD 

To remove the seasonal trends that are known to exist in the daily temperature data we 
need to determine the mean temperature for each day over all the years in the time series. 
We then calculate the fluctuation of the daily temperature from the mean daily temperature, 

AT, = Ti- <Ti> (1) 

where < > is the mean daily maximum temperature. Similarly we can also use the mean 
daily average temperature or the mean daily minimum temperature in place of the maximum 



temperature without changing the outcome of the analysis [28j]. The average temperature 
for some years can be higher or lower than the average temperature of the time series for 
the given location as a result of long-term atmospheric processes. To remove such remaining 
linear trends in the data we applied the Detrended Fluctuation Analysis (DFA) method , 
which removes trends in the data to allow for investigation of long-term correlations in the 
data. 

The noise and the nonstationarity in the temperature data usually hinders a reliable and 
direct calculation of the autocorrelation function C(n). Instead we calculate a running sum 
of the temperature fluctuations, 



y(m) = J2^T l (2) 



where m = 1, n. Next, the time series of the y(m) is divided into nonoverlapping intervals 
of equal length n. In each interval, we fit y(m) to a straight line, x{m) = km + d for each 



segment and calculate the detrended square variability F 2 (n) as 

(k+l)n 

F 2 (n)=<- (y(m) - x(m)) 2 > (3) 

m=fcn+l 

with 

N 

fc = 0,l,2,...,(--l). (4) 
If the temperature fluctuations were uncorrelated (white noise) we would expect 

F(n) ~ n a (5) 

where a = \. If a > |, we would expect long-range power law correlations in the data 
for the range of values considered. Fig. [1] shows one example of such an analysis method 
for New Bedford, MA, USA. The data set for this station spans 100 years of uninterrupted 
monthly temperature data between the years 1895 and 1994. The data follows a straight 
line with slope a = 0.66 ±0.01. Even though there is some scatter in the data after a period 
of 20 years, the scatter in the data is still within the error estimation of the analysis. This 
scatter is caused by the fact that the average in the detrended square variability F 2 (n) has 
been taken over larger and larger values of n, resulting in poor statistics. 

Fig. [2] shows a map of the power law exponents for the available 1184 stations in the 
continental US together with the locations of the stations. This map has been generated 
from the data using a running average method to generate a grid of 30 30 regions over 
the geographical region. A simple observation about the power law exponents is that even 
though the power law exponents crowd around the mean value of a = 0.60 ± 0.04, their 
spatial distribution is not uniform. From the map we can easily distinguish a few regions, 
namely the northeastern part is darker (lower power law exponent) than the southern and 
northwestern part (higher power law exponent). As in our previous work [43j, we cannot 
explain this behavior with either the distance from the oceans or the elevation or any simple 
dependence on geography. 

III. RESULTS 

To make some sense of this behavior we took a rather different approach. In the recorded 
human history there have been quite a few locally important climate changes. Many of 
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FIG. 1: The power law relationship between the time period and the detrended square variability 
for New Bedford, MA. 

these changes have been recorded in the paleoclimatic data. For example, tree ring history 
from Scandinavia indicates that there was a period of high temperatures between the 9th 
and 13th centuries, called the "Medieval Warm Period" 5(|. The tree ring temperature 
history for this period agrees with the glacier data, and is supported by the historical 
records of Norse seafaring and colonization around the North Atlantic at the end of the 
9th century. It is known that during this time the warmer climate helped in producing 



greater harvests in Iceland and parts of Greenland [51( which in turn helped the colonies 
in Greenland. But when we go to older times, it becomes more and more difficult to relate 
the paleoclimatic proxies to the climate of a certain region. From paleoclimatic proxies, like 
tree ring indices or stacked oxygen-isotope curves derived from deep sea cores, we normally 
obtain fluctuations of the local temperatures rather than their absolute values. Therefore, 
it is difficult to construct a method to characterize different climates based solely on the 
fluctuations of temperature values. As an example of this we have plotted in Fig. [3] the 
yearly average temperatures reconstructed from stalagmite thicknesses from Shihua Cave 



near Beijing, China 52j. Even though these average temperatures do not show any time 
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FIG. 2: The map of US with power law exponents and the stations marked. 

dependence, it is widely believed that the climate change in Central Asia around 100 BC 
was the main reason for the Hsiung-nu expansion towards the west |53|. Beijing, China is 
not considered a part of Central Asia, but it is a nearby location for which we have reliable 
data. When we plot the scaling exponents obtained over 400 year periods for this region 
together with the temperature data, we see that scaling exponents exhibit a jump around 
100 BC whereas the temperatures stay fairly constant (Fig. [3]). This leads us to believe 
that temperature fluctuations, together with the scaling exponents, can be used as climate 
indicators in the absence of other indicators like pan evaporation or precipitation. 



The Koeppen climate classification 54j is one of the most widely used climate classification 
systems. Although this system is based on the concept that native vegetation expresses the 
climate, it combines average annual and monthly temperatures and precipitation, and the 
seasonality of precipitation to classify climates. According to this classification scheme, 
there are about 14 different climate classes in the continental US. These classes range from 
tropical monsoon climate (Am) to hemiboreal climates (Dfb, Dsb, Dwb). Some of these 
classes are similar to each other in their characteristics but there is criticism because some 
regions are very broad. For example, Florida and New Jersey are both in the class Cfa 
(humid subtropical). In our work, our basic assumption is that the temperature fluctuations 
may be used to classify climates in the absence of other indicators. Therefore, we divided 
the climate of the US into distinct groups. These groups have tenuous connections to 
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FIG. 3: Temperature recontruction from the stalagmite thicknesses in Shihua Cave, Beijing, China 
(lines). Also plotted on this graph are the scaling exponents obtained from the same data for 
periods of 400 years (circles) . 

Koeppen classes. We then plotted the scaling exponents versus the standard deviation of 
the temperature fluctuations for locations fitting into these climate classes. Fig. H] shows 
these sample locations. As we can see from this graph, locations with similar climates tend 
to have similar scaling exponent / standard deviation values. 

We must be clear about one point: The reason why this method is being utilized is not 
to find an alternate classification scheme to already existing ones. The basic question to be 
answered is the following: When we have the temperature data from a paleoclimatic proxy, 
we know what the temperatures are, but we have little information about the climate. 
Can we make sense out of the temperature fluctuations from that location? Is it more 
like Ospina Perez, Colombia or closer to Las Vegas, NV? If we only have the standard 
deviations of the temperature fluctuations and the scaling exponents resulting from those 
distributions, obtaining clusters belonging to different climates would be extremely difficult. 
However, we possess additional information in the location of the weather stations. Therefore 
our task is not finding appropriate clusters, but classifying new data into already existing 
clusters. To aid us in classification and to support the claim that the differences in the 
scaling exponents of different geographical/climatic regions are statistically significant we 
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FIG. 4: The scaling exponents plotted against the standard deviation of the temperature fluctu- 
ations for different climate types. Examples: Black circles - Pensacola FL, black triangles - Ann 
Arbor MI, white squares - Minneapolis MN, Half-filled circles - Boulder CO, Grey triangles - San 
Antonio TX, Half-filled diamonds - Blacksburg VA, white circles - Stillwater OK 



have used the support vector machine (SVM) algorithm 55] for data classification. 

The inputs for a regular data classification task are training and testing data. The training 
dataset consists of one target value (in our case belonging to a specific climate type) and 
several features (like the standard deviation of the temperature fluctuations and the scaling 
exponent). We also need to tell the algorithm what we do not want. In the training set we 
have used 10 positive target values (belonging to the climate class we are analyzin g) a nd 20 
negative target points (belonging to different climate classes). The aim of SVM 56j is to 
produce a model which then predicts the target value of data instances in the test set. In 
the test set, we supplied the algorithm with all of the stations in our dataset and asked the 
program to identify the different climate types. The performance of such an algorithm is 
usually quantified by its accuracy during the test phase which mainly depends on the correct 
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treatment of true positives (TP) and true negatives (TN). It is also important to distinguish 
between two types of errors: A false positive (FP) and a false negative (FN). Consequently, 
the performance of the prediction is better judged if we add two more quantifiers, sensitivity 
and specificity. The accuracy of the data classification is defined as the ratio between the 
number of correctly identified samples and the total number of samples: 

TP + TN 

aCCUraCV = TP + FP + TN + FN (6) 

The sensitivity is the ratio between the number of true positive predictions and the 
number of positive instances in the test set: 

sensitivity = Tp — (7) 

Finally, the specificity is defined as the ratio between the number of true negative pre- 
dictions and the number of negative instances in the test set: 

TN 

specificity = Fp + TN (8) 

For our analysis, we used the different climates mentioned above in Fig. HI The results are 
summarized in Table [B In this table we see that our ability to predict accurately, sensitively 
and specifically never goes below 96%. However at this point we must explain why we came 
up with seven different climate classes. If we run the same classification algorithm with 
more than seven classes, the accuracy, or the sensitivity, or the specificity goes below 95 
%. Conversely, we can get even better accuracy/sensitivity/specificity by decreasing the 
number of classes. Therefore, we restricted ourselves to identifying seven different climatic 
regions. When dealing with paleoclimatic data, the accuracy/sensitivity/specificity levels 
can be adjusted to the needs of the situation. 

In Fig. [5] we show the result for all the available US data. Finally we show these locations 
on the US map (Fig. [6]). As can be expected, the result is not always perfect, however, we 
should keep in mind that these results have been obtained without other indicators such as 
precipitation or evaporation. The addition of such parameters into the same analysis system 
may make this tool a better categorization method, however as that area has been contested 
by many climatologists for more than a hundred years, we have restricted our analysis and 
results to the classification problem at hand. 
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TABLE I: SVM analysis of the 1184 weather stations in different climate zones. Accuracy, sensi- 
tivity, and specificity are defined in the text. 



Climate Type Accuracy Sensitivity Specificity 

Type 1 (Pensacola FL) 99.7% 100.0% 99.7% 

Type 2 (Ann Arbor MI) 98.1% 96.0% 99.3% 

Type 3 (Minneapolis MN) 98.8% 99.3% 98.7% 

Type 4 (Boulder CO) 99.3% 100.0% 99.3% 

Type 5 (San Antonio TX) 97.8% 99.3% 97.6% 

Type 6 (Blacksburg VA) 98.6% 100.0% 98.4% 

Type 7 (Stillwater OK) 96.0% 96.4% 96.0% 



IV. CONCLUSION 



We have investigated the power law behavior of the temperature fluctuations to gain 
insight into a statistical description of different climate types. Normally we need other 
indicators such as the pan evaporation and precipitation to classify the climates, but two 
different fields, both connected with the global climate change, led us to this approach. 
First, different projections for the future of the climate change give us valuable data about 
the future temperatures and possibly precipitation. And second, paleoclimatology gives us 
different proxies to determine the temperatures of ancient times. However, as both of these 
areas are focused on determining the temperatures, our aim is to aid this area of research 
by developing a tool that uses temperature fluctuations to determine the climate. 

Our results suggest that on the basis of the power law scaling of the temperature fluc- 
tuations we can distinguish between different climates. This is by no means a new way to 
classify climates. It is only intended to compare the climate of a spatio-temporal location 
for which we only have temperature fluctuations without the exact units to the present day 
climate for which we possess other meteorological measurements in addition to temperature. 
The real challenge comes when we try to expand this analysis into the past, since, to our 
knowledge, no reliable monthly data exists beyond 218 years ago 23[]. However, previous 
work where authors have established long-range power law behavior using rescaled-range 
analysis [57], [58( gives us hope about expanding the use of this method we have developed. 
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FIG. 5: The scaling exponents plotted against the standard deviation of the temperature fluctu- 
ations for different climate types after the classification by the SVM algorithm. Examples: Black 
circles - Type 1, black triangles - Type 2, white squares - Type 3, Half-filled circles - Type 4, Grey 
triangles - Type 5, Half-filled diamonds - Type 6, white circles - Type 7 

The climate projections for the future are crucial in our understanding of how climate 
change impacts the environment. When we use these global climate models to extrapolate 
the climate changes into the future, the method of Detrended Fluctuation Analysis proves to 
be a powerful tool 



19 



30j. The present work may help to fill in the voids in existing climate 



models by providing further information on climate and temperature while improvements 
are made on the length scales these climate models cover along with their applicability to 
future timescales. 
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FIG. 6: The map of US with the classification results for different climate types and the stations 
marked. Examples: Black circles - Type 1, black triangles - Type 2, white squares - Type 3, 
Half-filled circles - Type 4, Grey triangles - Type 5, Half-filled diamonds - Type 6, white circles - 
Type 7 

Acknowledgments 

Data was provided by U.S. Historical Climatology Network and the National Climatic 
Data Center. 



[1] Global Glacier Recession, National Snow and Ice Data Center (2006), 

http:/ / nsidc.org/ glims / glaciermelt / index.html. 
[2] M. S. Pelto, The Cryosphere Discuss. 1, 17 (2007). 

[3] Antarctic ice shelf breaks apart, BBC News (2002), http://news.bbc. co.uk/2/ 

hi /science /nature/1880566. stm. 
[4] Ice shelf in warmest part of Antarctica collapses, USA Today (2003), 

http:/ / www.usatoday.com / news / science / cold-science/2003-03- 19-ant-ice.htm. 
[5] Larsen B Ice Shelf Collapses in Antarctica, National Snow and Ice Data Center (2002), 

http:/ /nsidc.org/iceshelves/larsenb2002/. 



13 



[6] T. L. Root, J. T. Price, K. R. Hall, S. H. Schneider, C. Rosenzweig, and A. Pounds, Nature 
421, 6918 (2003). 

[7] T. L. Root, D. MacMynowski, M. D. Mastrandrea, and S. H. Schneider, PNAS 102, 7465 
(2005). 

[8] C. D. Thomas, A. Cameron, R. E. Green, M. Bakkenes, L. J. Beaumont, Y. C. Collingham, 

B. F. N. Erasmus, M. F. de Siqueira, A. Grainger, L. Hannah, et al., Nature 427, 145 (2004). 
[9] H. C. Fritts, Tree Rings and Climate (Academic Press, 1976). 
[10] J. R. Petit, J. Jouzel, D. Raynaud, N. I. Barkov, J. M. Barnola, I. Basile, M. Bender, J. Chap- 

pellaz, J. Davis, G. Delaygue, et al., Nature 399, 429 (1999). 
[11] T. M. Hill, J. P. Kennett, D. K. Pak, R. Behl, C. Robert, and L. Beaufort, Quaternary Science 

Reviews 25, 2835 (2006). 
[12] F. Gingele, P. de Deckker, and M. Norman, Earth and Planetary Sci. Lett. 255, 257 (2007). 
[13] R. B. Dunbar, G. Wellington, M. Colgan, and P. Glynn, Paleooceonography 9, 291 (1994). 
[14] J. W. Kantelhardt, D. Rybski, S. Zschicgner, P. Braun, E. Koscielny-Bunde, V. Livina, 

S. Havlin, and A. Bunde, Physica A 330, 240 (2003). 
[15] J. Eichner, E. Koscielny-Bunde, A. Bunde, S. Havlin, and H. J. Schnellnhuber, Phys. Rev. E 

68, 046133 (2003). 

[16] V. Livina, Y. Ashkenazy, P. Braun, R. A. Monetti, A. Bunde, and S. Havlin, Phys. Rev. E 

67, 042101 (2003). 
[17] R. A. Monetti, S. Havlin, and A. Bunde, Physica A 320, 581 (2003). 
[18] A. Bunde and S. Havlin, Physica A 314, 15 (2002). 

[19] R. Govindan, D. Vjushin, A. Bunde, S. Brenner, Y. Ashkenazy, S. Havlin, and H. J. 
Schnellnhuber, Phys. Rev. Lett. 89, 028501 (2002). 

[20] D. Vjushin, R. Govindan, S. Brenner, A. Bunde, S. Brenner, S. Havlin, and H. J. Schnellnhu- 
ber, J. Phys.: Cond. Matt. 14, 2275 (2002). 

[21] J. W. Kantelhardt, E. Koscielny-Bunde, H. H. A. Rego, S. Havlin, and A. Bunde, Physica A 
295, 441 (2001). 

[22] R. Govindan, D. Vjushin, S. Brenner, A. Bunde, S. Havlin, and H. J. Schnellnhuber, Physica 
A 294, 239 (2001). 

[23] E. Koscielny-Bunde, A. Bunde, S. Havlin, H. E. Roman, Y. Goldreich, and H. J. Schnellnhuber, 
Phys. Rev. Lett. 81, 729 (1998). 

14 



[24] E. Koscielny-Bunde, A. Bunde, S. Havlin, and Y. Goldreich, Physica A 231, 231 (1996). 

[25] L. Bodri, Theoretical and Applied Climatology 51, 51 (1995). 

[26] L. Bodri, Theoretical and Applied Climatology 49, 53 (1994). 

[27] R. O. Weber and P. Talkner, J. Geophys. Res. Atmos. 106, 20131 (2001). 

[28] P. Talkner and R. O. Weber, Phys. Rev. E 62, 150 (2000). 

[29] R. O. Weber, P. Talkner, and G. Stefanicki, Geophys. Res. Lett. 21, 673 (1994). 
[30] K. Fraedrich and R. Blender, Phys. Rev. Lett. 90, 108501 (2003). 
[31] K. Ivanova and T. Ackermann, Phys. Rev. E 59, 2778 (1999). 
[32] G. Kiely and K. Ivanova, Phys. Chem. Earth B 24, 781 (1999). 
[33] K. Ivanova and M. Ausloos, Physica A 274, 349 (1999). 

[34] K. Ivanova, M. Ausloos, E. Clothiaux, and T. Ackermann, Europhys. Lett. 52, 40 (2000). 
[35] M. Ausloos and K. Ivanova, Phys. Rev. E 6304, 047201 (2001). 

[36] K. Ivanova, E. Clothiaux, H. Shirer, T. Ackermann, J. Liljegren, and M. Ausloos, J. Appl. 

Meteorology 41, 56 (2002). 
[37] N. Kitova, K. Ivanova, M. Ausloos, T. Ackermann, E. Clothiaux, and M. Mikhalev, Int. J. 

Mod. Phys. C 13, 217 (2002). 
[38] K. Ivanova, H. Shirer, E. Clothiaux, N. Kitova, M. Mikhalev, T. Ackermann, and M. Ausloos, 

Physica A 308, 518 (2002). 
[39] K. Ivanova and M. Ausloos, J. Geophys. Res. - Atm. 107, 4708 (2002). 

[40] K. Ivanova, T. Ackermann, E. Clothiaux, P. Ivanov, H. Stanley, and M. Ausloos, J. Geophys. 

Res. - Atm. 108, 4268 (2003). 
[41] M. Ausloos, Physica A 336, 93 (2004). 

[42] M. L. Kurnaz, J. Stat. Mech.: Theory and Experiment p. P07009 (2004). 
[43] M. L. Kurnaz, Fractals p. 365 (2004). 

[44] B. B. Mandelbrot and J. R. Wallis, Water Resources Research 5, 321 (1969). 

[45] B. B. Mandelbrot and J. R. Wallis, Water Resources Research 5, 967 (1969). 

[46] A. Arneodo, E. Bacry, P. V. Graves, and J. F. Muzy, Phys. Rev. Lett. 74, 3293 (1995). 

[47] C. K. Peng, S. V. Buldyrev, S. Havlin, M. Simmons, H. E. Stanley, and A. L. Goldberger, 

Phys. Rev. E 49, 1685 (1994). 
[48] K. Fraedrich and R. Blender, Phys. Rev. Lett. 92, 039802 (2004). 
[49] R. Blender and K. Fraedrich, Geophys. Res. Lett. 30, 1769 (2003). 

15 



[50] R. Cook, Climate Dynamics 11, 211 (1995). 

[51] E. R. Bryson and T. J. Murray, Climates of Hunger (University of Wisconsin Press, Madison, 
1977). 

[52] M. Tan, T. S. Liu, J. Hou, X. Qin, H. Zhang, and T. Li, Geophysical Research Letters 30, 
1617 (2003). 

[53] C. Chase-Dunn, R. Niemeyer, A. Alvarez, H. Inoue, K. Lawrence, and A. Carlson (2006). 
[54] W. Koeppen, Geographischen Zeitschrift 9, 593 (1900). 

[55] V. Vapnik, The Nature of Statistical Learning Theory (Springer Verlag, 1995). 

[56] T. Joachims, Advances in Kernel Methods - Support Vector Learning (M.I.T. Press, 1999). 

[57] R. H. FluegemanJr and R. S. Snow, Pageoph. 131, 307 (1989). 

[58] Y. Wang, M. E. Evans, T. C. Xu, N. Rutter, Z. Ding, and X. M. Liu, Can. J. Earth Sci. 29, 
296 (1992). 



16 



